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ABSTRACT 
A  probabilistic  model  is  developed  for  treating  the 
problem  of  optimizing  the  intended  trajectory  and  associated 
height  of  burst  of  a  missile-warhead.   The  angle  of  re-entry 
is  treated  as  a  control  variable.   Computational  techniques 
which  may  be  operationally  acceptable  are  described. 
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I   INTRODUCTION 

Consider  the  problem  of  setting  the  fusing  mechanism  of  a 
warhead  aboard  a  ballistic  missile  in  order  to  optimize  the 
probability  of  destroying  a  target. 

The  general  problem  as  stated  above  becomes  extremely 
complicated  when  one  considers  such  things  as  missile  in-flight 
failures,  warhead  failures,  back-up  fusing  mechanisms  employed, 
type  of  warhead  and  type  of  burst,  distance  from  launcher  to 
target  (which  will  affect  range  and  Azimuth  errors),  geodesic 
error,  (i.e.,  to  what  degree  of  accuracy  do  we  know  the 
distance  to  target.) ,  and  more  frustrating,  the  problem  of 
partial  kills,  i.e.,  the  partial  destruction  of  the  target  or 
a  portion  of  the  target  for  a  period  of  time  less  than  the 
duration  of  hostilities. 

In  order  to  reduce  the  problem  to  something  that  can  be 
handled,  it  will  be  convenient  to  cover  briefly  three  topics: 
the  concept  of  ballistic  coefficient,  the  reflected  (mach) 
wave  phenomenon  and  the  definite  range  law  (cookie  cutter) . 
Ballistic  Coefficient: 

For  our  purposes,  we  may  think  of  the  ballistic  coef- 
ficient as  a  parameter  which  describes  how  much  a  re-entry 
vehicle  will  be  affected  by  the  atmosphere  (see  Fig.  1).   It 
is  proportional  to  the  ratio  of  the  weight  of  the  vehicle  to 
its  area,  i.e.,  the  square  of  the  maximum  diameter  (see  Ref.  4); 
e.g.,  for  a  given  area,  an  increase  in  weight  will  tend  to 
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Fig.  I  —  Trajectory  contours  for  10-deg  re-entry  angle 


"straighten"  the  trajectory  whereas  a  decrease  in  weight  will 
result  in  a  slower  re-entry  and  thereby  increase  the  time  spent 
in  the  gravitational  field.   The  result  will  be  a  more  curved 
trajectory. 
Mach  Effect: 

When  a  warhead  is  burst  in  the  air  close  to  ground  (e.g., 
300  feet  for  a  1  KT  warhead),  there  is  a  region  of  unusually 
high  overpressures  due  to  the  merging  of  the  incident  and  re- 
flected waves  (see  Ref.  1  and  Fig  2).   This  region  is  known  as 
the  Mach  Region.   Because  of  this  the  curves  in  Fig.  2  are  not 
monotonic  but  have  the  unusual  shape  exhibited  there. 
Definite  Range  Law: 

This  is  the  "cookie-cutter"  concept  that  is  frequently 
used.   Because  of  its  simplicity  it  seems  to  give  useful 
results  for  many  objectives. 

A  couple  of  examples  should  be  sufficient  for  our  purposes. 
Consider  a  point  target  in  space  at  which  we  are  firing  a  pro- 
jectile with  a  kill  range  of  one  mile.   Imagine  the  target  to 
be  at  the  center  of  a  sphere  of  radius  one  mile.   If  we  detonate 
this  weapon  anywhere  within  the  sphere  we  consider  the  target 
destroyed;  anywhere  outside  the  sphere  is  considered  a  miss.   Of 
course,  no  weapon  behaves  in  such  a  manner.   "Partial  kills"  are 
generally  the  rule  rather  than  the  exception,  and  they  are  very 
difficult  to  quantify.   The  definite  range  law  is  a  mathematical 
convenience.   Presumably,  this  range  can  be  regarded  so  that  a 
controlled  fraction  of  partial  kills  will  be  treated  as  misses 
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T^g.- 2~ repYesents  a  destruction -curve  for  various  overpr'es'su^es/  .of  a '  ,  '  •  ^7' 
1  KT  weapon.  For  a  weapon  of  different;'  yield,  W  ^the'/ording te,  and  'abscissa  '  ' i 
are  each  multiplied  by  W-*.'.  See  Ref.  1.  For  example^'/a  '600.'  foot  high  burst 
of  a  1  KT  weapon,  the  20  P.S.I,  curve  will',  ex  texid  to  about  940 1  feet' from  :  ' 
ground  zero.  Similarly,  a  1  MT  weapon  det'o'natedi  at  .'O.OOO'/feet  will  have  a  ;,' 
20  P.S.I,  overpressure  to  about  9,400.  feet.  ■   ''  '   ,v  ( 


(and  a  controlled  fraction  will  be  scored  as  hits). 

As  another  example  consider  Fig.  5.   If  we  imagine  a  target 
at  (0,0)  and  estimate  that  it  takes  about  20  P.S.I,  of  over- 
pressure for  destruction,  the  20  P.S.I,  curve  shown  could  be 
considered  as  a  "cooker-cutter";  i.e.,  a  burst  just  inside 
will  destroy  the  target  and  one  just  outside  will  not. 

In  light  of  the  above  discussion,  the  following  three 
simplifying  assumptions  will  serve  as  the  basis  of  this  paper: 

1.  The  re-entry  vehicle  has  a  sufficiently  high  ballistic 
coefficient  so  that  the  trajectory  may  be  assumed  to  be  a 
straight  line  in  the  neighborhood  of  the  target.   Moreover, 
the  vehicle  is  launched  from  a  sufficiently  great  distance  so 
that  all  possible  paths  may  be  assumed  to  be  parallel  in  the 
neighborhood  of  the  target. 

2.  The  overpressure  contours  of  Ref.  1  (and  hence  of  Fig.  2) 
are  valid  and  the  target  damage  is  a  function  of  overpressure 
alone . 

3.  A  definite  range  law  holds  and  20  P.S.I,  is  necessary  for 
destruction. 

The  general  objective  is  to  select  an  optimum  height  of 
burst.   Further  examination  of  the  problem  reveals  that  some 
other  parameters  enter  into  the  picture.   The  angle  that  the 
vehicle's  trajectory  makes  with  the  surface  can  be  selected 
in  advance.   There  is  an  advantage  to  small  angles  since  then 
the  missiles  are  more  difficult  to  detect  (and  hence  defend 
against).   Larger  angles  generally  lead  to  higher  kill 
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probabilities.   The  capability  to  vary  this  parameter  will  make 
defense  more  difficult.   Given  a  trajectory  angle,  the  optimum 
burst  point  becomes  a  function  of  height  of  burst  and  aim  point 
(intersection  of  the  intended  trajectory  with  the  surface) . 
Locating  the  target  at  (0,0,0)  and  using  the  Cartesian  system 
wi  th 

)(,=  downrange  distance 
Yx=   crossrange  distance 
A3=  vertical 
and  introducing  parameters 

-0-=   angle  the  trajectory  makes  with  the  surface 
S  =  distance  from  aim  point  to  target 
M  =  intended  height  of  burst 
(see  Fig.  3),  the  problem  now  may  be  stated: 

Given-^  ,  choose  the  pair  (  5  ,  LK.  )  so  that  the  probability 
of  target  destruction  is  maximized. 

This  paper  must  be  viewed  as  a  pilot  study  for  solving  the 
problem  outlined  above.   It  will  be  noted  that  a  large  amount  of 
detailed  work  remains  to  be  done  in  order  to  provide  good  answers 
The  methods  herein  can  serve  as  outlines  and  guides. 

One  of  the  more  important  results  is  the  development  of  a 
graphical  computation  system  using  dividers,  graphs  of  over- 
pressure functions  and  correction  curves.   Thus,  the  optimum 
solution  may  be  obtained  under  operational  conditions  without 
the  need  of  a  digital  computer. 


The  organization  of  the  report  is  as  follows:   The  mathe- 
matical model  and  nature  of  the  problems  of  mathematical 
analysis  are  presented  in  section  II.   The  computational  problems 
are  described  in  section  III.   The  operational  hand  approxi- 
mations and  corrective  curves  are  given  in  section  IV.   The 
conclusions  appear  in  section  V.   The  Monte  Carlo  technique 
employed,  along  with  the  computer  program,  appear  in  Appendices 
A  and  B  respectively.   Appendix  C  contains  a  tabulation  of 
the  parameter  values  determined  by  the  graphical  method  de- 
veloped in  section  IV. 


II 


CONSTRUCTION  OF  THE  PROBABILISTIC  MODEL; 
MATHEMATICAL  CHARACTER  OF  THE  PROBLEM 


We  have  already  introduced  the  coordinate  system;  the 
angle  the  trajectory  makes  with  the  surface;  the  intended 
height  of  burstM   ;  the  aim  point  (5  ,  0,  0).   The  trajectory 
is  actually  a  random  phenomenon  and  by  assumption  1  can  be 
characterized  by  the  impact  point  and  the  incidence  angle  -Q~   . 
The  actual  height  of  burst  is  also  subject  to  random  errors. 
Thus,  we  let 

(Y^ ,  Y2 ,  0   )  be  the  random  vector  representing  the  impact 
point , 

(X\ ,   X2 ,  X3  )  be  the  random  vector  representing  the  burst 
point 
A  sample  situation  is  given  in  Fig.  4. 

It  is  assumed  that: 
!•   Y]_,Y2has  a  circular  normal  distribution  with  mean  vector 
£  ,  0   and  covariance  matrix  0\,  I  where  I  is  the  identity 

matrix. 

2 
2.   X3  has  a  normal  distribution  with  mean U  ,    variance^  and  is 

independent  of  Yi ,  and  \. 

From  Fig.  4  it  is  seen  that  Xi  =  Yi  -  X3  cot  "6" 

X2  =  Y2 


It  follows  that  the  probability  density  function  of  the 

detonation  point  is  ,  ~  -y-, 

fetfT*<rcp  v 

The  ballistic  error  parameter  07  ,  and  the  fusing  error  parameter 
(T^'should  be  determined  from  independent  experiments.   There  are 
many  types  of  warheads  and  fuses  and  we  merely  carry  (P  ,  and  (ff 
as  "known"  parameters. 

Since  20  P.S.I,  will  destroy  the  target,  the  event  of 
target  destruction  can  be  visualized  in  the  (Xx ,  X2,  X3)  space 
by  rotating  the  curve  in  Fig.  2  about  the  vertical  axis.   If 
we  call  this  set  A,  then  the  probability  of  target  destruction 
may  be  expressed  as  the  triple  integral  of  the  density  function 
(2.1)  over  the  set  A 

(2,2)      Pd  =    fffffcs&.XsjJr.AJrj 

It  may  be  well  to  clarify  a  small  logical  difficulty.   The 
probabilistic  model  permits  negative  values  for  the  height  of 
burst.   These  are  physically  impossible  and  should  be  interpreted 
as  a  duded  warhead  that  impacts  the  surface  without  exploding. 

The  optimum  solution  can  be  obtained  by  taking  the  appro- 
priate partial  derivatives  of  (2.2),  equating  them  to  zero  and 
solving  the  resulting  system  of  equations.   Thus 

dp 


(2.3) 


SJH  *  ffffx,  -ju)f(W>)**  fc  4*j 


77  *  fffftrf+XiCote-ffa^jJxM 

A        10 


Examination  of  the  character  of  equations  (2.4)  is  facili- 
tated if  we  make  the  change  of  variables 

(2.5)    y,    -   X U.   ,  yz-  _   ,$-->£ 

and  let 

In  this  coordinate  system  the  set  A  is  transformed  into  a  set  B 
the  exact  nature  of  which  is  difficult  to  visualize  (see  Fig.  2), 
but  we  can  say  B  is  bounded  since  the  transformation  is  linear. 
In  the  new  system,  equations  (2.4)  may  be  expressed  as 

Geometrically,  the  problem  now  may  be  viewed  as  follows:   Given 
a  set  B,  choose  a  centering  point  (  J  ,'-'•' ,/'0  so  that  the 
probability  of  B  is  maximized. 

All  solutions  of  (2.6)  are  critical  points  which  may  be 
local  maxima,  minimum  or  saddle  points.  The  nature  of  these 
points  may  be  examined  by  means  of  second  order  derivatives. 

LettinE   ,,  ,,*  -tffvi-dV.Vtf"?? 

we   can  write  i 

(2.8)     ^--JJTtOS-^-^^^.^^'^^ 

(2.9)  p/Mf  =  f(f  [(Wi'tyM*.  \,ii)  H**A* 


li 


A  critical  point  will  be  a  local  maximum  if 

(2-n)  %a  %    '  %  ><*      and  da^<  ^        (Ref-  3  p-  232) 

It  can  generally  be  said  that  Daam  (and  D,  -   ,  for  the  same 


nee  ffr(%-s>y(iw)Wii! 


reason)  is  negative  since  the  variai 
of  Y3  is  one  and  (2.8)  compares  a  truncated  (by  B)  version  of 
this  variance  with  unity.   The  remaining  part  of  (2.11)  is  a 
difficult  question,  however,  and  may  require  more  specific  in- 
formation about  the  set  B. 

The  set  B  plays  a  role  in  determining  the  number  of  critical 
points.   For  example,  if  B  were  a  sufficiently  long  dumbell 
shaped  region,  there  would  be  at  least  two  solutions  to  (2.6). 
It  would  be  interesting  to  examine  the  question  of  number  of 
critical  points  if  B  were  convex.  With  this  condition  it  may  be 
possible  to  show  there  is  only  one.   This  contingency  was  not 
examined  since  the  set  A  (and  hence  B)  is  not  convex. 

Since  the  set  A  is  formed  by  rotation  of  Fig.  2  about  the 
X3  axis,  it  can  be  shown  that 

(2. 12)     ttfxx-f  (*, t  ^  *, )  dx{  dxx  J/j  -  0 

n 

The  same  must  be  true  in  the  transformed  coordinate  system,  and 

it  follows  from  (2.6)  and  (2.12)  that  (J  ,  0,iU)  will  be  the 
center  of  gravity  of  the  conditional  probability  density  given 
(Y]_,Y2  ,Y3  )  belong  to  B  (if  J  and  M  are  unique  solutions  of 
(2.6)). 
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Ill 

COMPUTATIONAL  PROBLEMS 

Computer  techniques  are  required  for  the  solution  of  the 
system  of  equations  (2.4) ,  the  examination  of  the  critical 
points  (2.11),  and  the  determination  of  the  optimum  probability 
of  target  destruction  (2.2).   It  was  decided  to  use  Monte  Carlo 
techniques  (see  Appendix  A)  in  this  study,  and  the  first  step  is 
to  characterize  the  set  A.   The  curve  in  Fig.  2  may  be  viewed 
as  giving  distance  as  a  function. 

This  curve  is  not  given  in  analytic  form.   Also,  note  of 
height,  i.e.,  g(Xo  )  that  there  is  a  small  interval  of  X3 
values  for  which  g(Xo  )  is  double  valued.   It  was  decided  to 
ignore  the  double  valued  feature  and  fit  a  polynomial   of 
order  7  by  the  method  of  least  squares,  using  23  values  read 
directly  from  the  graph.   This  fit  appears  in  Fig.  8.   Using 

the  polynomial  as  g(X3  ),  the  set  A  may  be  represented  by 

f  2    2      2 

(3.1)   A  =  |  X1,X2,X3  :   Xx+  X2  £  g  (X3)l 

The  integrals  in  (2.2),  (2.4)  and  (2.11)  all  are  similar 
in  nature.   A  rapidly  converging  iterative  technique  may 

generally  be  available  for  the  solution  of  the  system  (2.3). 

LYX 
Letting     y  ,    -  2  * 

the  equations  (2.4)  may  be  written 

°  =  tffarfUt  cor  6-)  p(te±^*)0(£)fl&Jd)rMJ* 

A  13  ' 


Let  functions  gl(/U/  J  )  and  g2(/vv ,  /  )  be  defined  so  that  the  above 
equations  are 

Then,  using  an  initial  approximation  of  AA0      f        we  can  use 
the  iterative  scheme 

(3-3)   /"*♦,--?,  (^,,  /J  ,  fH,*?*fa?*) 

provided  it  converges. 

The  examination  of  the  question  of  convergence  usually 
goes  along  the  following  lines:   If  the  series 


(3.4) 

converges  absolutely,  then  AA^  and  f    converge  to  finite 
numbers  which  will  be  the  solutions  of  2.4.   This  will  take 
place  if  a  constant  r  (  0<r<l)  can  be  found  such  that 

for  then  - 

The  question  of  convergence  of  the  scheme  is  more  easily 
examined  after  making  the  transformation  (2.5).   Let  this  be 
done  and  write   P*  -  tf {  0(ir  f) /$  W  H  U,  ~/*)  <**>  ^  *  ^  % 
and  then  . 

(,.6)    ?,fa.sj--£U[yif(yl-f*)j<vftw> 

o 
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Proceeding  as  indicated, 
(3.7) 


From  the  theorem  of  the  mean 

0(tl-  0(f)  i  U-t'l  ^wf(f) 

it  is  easily  shown  that 

and  this  value  may  be  used  in  3.7.   Thus  we  can  write 

|  4  (/, -/*,)- 4  (/, -/*,.,)!  *  P \h /V-< I 
\0(Y,-f,)-^Y.-f/.,)lt(>\?l-  ?ril 

Sharper  bounds  may  be  available  if  the  maximum  is  3.8  and  the 
points  (  F  ,0  >/^«)  are  constrained  to  be  in  the  set  B.   The 
second  term  of  3.7  can  be  estimated  in  a  similar  fashion: 

which  can  be  treated  as  above. 
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Collecting,  we  have  , .  ,   , 

(3.10) 

which  will  have  the  form  (3.5)  if  it  can  be  shown  that  the 
coefficients  are  less  than  one.   This  is  a  difficult  analytic 
problem  and  is  beyond  the  scope  of  this  paper.   The  technique 
(3.3)  was  used  in  our  computations  and  lend  to  answers  that 
were  stable  in  the  light  of  our  capability  to  compute  by  Monte 
Carlo.  The  signs  of  the  determinants  (2.11)  were  computed  at 
the  solution  values  and  all  were  found  to  be  positive. 
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ACTUAL   BURST 


FIGURE  4 
SKETCH  OF  TRAJECTORY  AND  CO-ORDINATE  SYSTEM 
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IV 

A  SIMPLIFIED  APPROACH  TO  OBTAIN  AN  APPROXIMATE  SOLUTION; 
GRAPHICAL  METHODS  AND  CORRECTION 


If  the  down  range  and  cross  range  errors  are  assumed  to 
be  small,  the  problem  can  be  reduced  to  a  simple  two-dimen- 
sional model  since  we  intuitively  feel  that  a  'best'  trajectory 
is  one  that  passes  directly  over  target. 

At  this  point  we  will  assume  we  follow  our  intended 
trajectory  with  probability  1  and  find  values  of  AA   and  >  that 
optimize  this  conditional  probability  of  destruction. 

If  Fig.  2  were  reflected  about  a  vertical  axis,  Fig.  5 
will  result.   Again,  assume  the  target  to  be  at  (0,0)  and  we 
wish  to  detonate  our  weapon  inside  the  20  P.S.I,  envelope. 

We  will  make  the  following  assumption  to  supplement  those 
in  section  III: 

The  warhead  passes  directly  over  target  in  a  straight 
line  trajectory  making  an  incidence  angle  of  0"with  the  surface. 
The  collection  of  possible  trajectories  will  intersect  the  20  P.S.I 
envelope  as  shown  in  Fig.  5.   Denote  the  lower  intercept(of  a 
given  trajectory)  by  a  and  the  upper  by  ©(  . 

The  maximization  takes  place  in  two  stages.   First,  holding 
■0  and  J  fixed,  find  the  best  value  of  A/(  as  a  function  of  J   . 
Second,  vary  )  until  the  maximum  is  obtained.   To  implement  this 
procedure ,  let  .     ,  ■> 


(4.1)   f  (M)   r  — '—       f    n     l<rl  Jv 
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Taking  the  derivative  and  equating  to  zero 


Making  the  substitution  t  =  Y-u 

r  ±£  £ 

j^  iff  T?      Jf-o 

<r 


*-  -0 


<r 


Since  a  function  and  its  log  have  a  miximum  at  the  same  points 


(*29 '-  (¥) 


Simplifying  and  solving  for  u 


I    ~  "<- 


/A  -    £s3 
(4-3)  /*«**=     —  ^?C° 

The  problem  is  degenerate  for  the  case  o(  =  a    since  the 
probability  of  detonation  is  zero.   It  is  interesting  to  note 
that  the  above  solution  is  independent  of (T  and  the  velocity 
of  re-entry. 

To  verify  we  have  a  maximum,  consider 
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Making  the  substitution  t  =  Y  -  u  and  simplifying  we  obtain 
after  one  integration  by  parts  we  see 


-  <**>  - (^  c 


which  is  true  for  all  *(  }  Q* 

The  value  of  F  at  the  point  Aimax   which  is  the  probability 
of  destruction  given  a  particular  line  is  given  by 


nwZ*r  €*»■■*> 


Le 


tting  ^  _         y  -    «^£fl 


FfW  = 


> 


,/rr      tC  Oi  A 


1(T 


*    ^5f 


or 


.x1 


(*•»  FW*  M  ^fee^f-i.o 


..  ^ 


So  far  we  have  the  optimum  height  of  the  IBP,  but  still 
must  find  out  where  to  plan  the  detonation,  i.e.,  to  select  f  . 
Since  we  are  interested  in  maximizing  the  probability  of 
destruction  (4.5),  we  see  this  is  clearly  equivalent  to 
maximizing  (  d^  -a)  since  the  function  is  monotonic.   Maximizing 
(  o(  -a)  is  in  turn  equivalent  to  finding  the  longest  line  that 
can  be  drawn  through  Fig.  5  at  an  angle  &  . 

In  Fig.  5,  a  series  of  parallel  lines  were  drawn  roughly 
every  20  feet  at  an  angle-©-  to  the  horizontal  and  the  longest 
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line  determined  by  the  use  of  dividers,  and  choosing  the  one 

C  7  ^ 

with  max   y  o^  -  C(    [      ,  If  a  maximum  occurred  between  lines  it 

was  deemed  to  be  at  the  mid-point  of  the  two  which  gives  an 

"accuracy"  of  10  feet  for  a  1  KT  weapon  or  100  feet  for  a  1  MT 

weapon . 

By  assuming  the  measurements  were  accurate  in  the  sense 

that  one  can  differentiate  between  a  line  being  longer  than  an 

adjacent  line  (which  is  reasonable  since  the  vast  majority 

were  easily  discernable)  we  can  say  with  probability  one,  the 

'/J 

maximum  error  is  10  W    feet  for  a  W  KT  weapon. 

Measurements  were  taken  in  this  manner  for  several  values 
of -O"  and  the  resulting  values  of  ilmax  were  plotted  against  S~    • 
Readings  were  concentrated  at  points  where  the  curve  appeared 
to  be  rapidly  changing.   A  similar  procedure  was  used  for  values 
of  ©■  vs  /  and  0-   vs  f     .   in  order  to  interpolate  polynomials 

of  degree  1,  2,  ,  20  were  fit  in  a  least-squares  sense  to 

the  data  points.   The  degree  chosen  from  the  one  with  the 
smallest  \    /  V  -  Q    I    >  where   / ^   =  emperical  value 

and  Oi  =  computed  value.   The  above  criterion  was  chosen  since 

it  was  planned  to  interpolate  values  of  •0-=5,  6,  7,  ,  90 

and  the  measurements  were  felt  to  be  fairly  accurate  to  start 
with  and  large  deviations  in  the  fit  would  only  tend  to 
aggrevate  the  original  data  points. 

The  polynomials  were  fed  into  the  computer  and  evaluated 
at  the  points  "©"=  5,  6,  90.   The  results  were  tabulated 

22 


in  table  1  where: 

Jimax  =  planned  burst  point  measured  vertically  along  HOB  axis 
XCOORD  =  planned  burst  point  measured  horizontally  along 
"DISTANCE  FROM  GROUND  ZERO"  AXIS 
-0*=  Angle  trajectory  intersects  ground  measured  from  the 

horizontal 
J  =  Aim  point  measured  positively  to  the  right  along  "DISTANCE 
FROM  GROUND  ZERO"  AXIS 

REMARKS  CONCERNING  THE  GRAPHICAL  APPROXIMATION 

The  most  interesting  discovery  is  that  the  'best'  IBP  is 
never  over  target;  for  low  angles  of  re-entry,  the  burst  is 
planned  prior  to  target;  for  steep  angles  of  re-entry,  the  IBP 
is  past  target. 

The  sharp  discontinuity  in  XCOORD  at  about  28°  is  caused  by 
entry  into  the  envelope  in  the  region  (850,  700)  of  Fig.  5  where 
the  slope  of  the  envelope  and  the  slope  of  the  trajectory  are 
nearly  equal . 

We  will  next  apply  the  correction  discussed  in  section  III. 
The  method  developed  above  was  used  as  the  first  approximation 
for  the  iterative  technique  developed  in  section  III.   The 
details  follow. 
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L«  P  -   J  i  J    l*it)*C<r  X. 

r-  ff 

Solving  for  AA  and   «: 

let    >CH ,  ,     F ,  =   First  approximation  obtained  by  using  the 

method  of  section  IV 

In  general,  about  10  iteratives  were  needed  f  or  ©>  >   10° 
and  about  20  for#  =  10°. 

Data  was  collected  for  -©•  =  10°,  20°,  -----80°  for  two 
sets  of  the  standard  deviation  (  (H  (JT  )  .   Table  1  contains  the 
percent  of  error  in  M  and  f  and  Fig.    contains  a  plot  of $-  vs 
percent  of  error. 

Thus,  operational  procedures  can  be  developed  that  do  not 
require  a  high  speed  computer.   The  first  approximation  is 
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<y  =  500,      <r\  m  200  6"  =  200,^  =  500 

e-  h  T  /*  f 


10 

+20.8 

+0.1 

+50 

+25 

20 

+  7.5 

+1.34 

+48.9 

+71 

30 

+  4.28 

-1.8 

+12.5 

+  5.9 

40 

+  7.9 

-  .134 

+  3.54 

-24.2 

50 

+  1.5 

-4.8 

-  6.15 

-45 

60 

+   .065 

-2.05 

-  5.6 

-41.5 

70 

+   .43 

-1.9 

-  7.1 

-65 

80 

+  1.3 

-3.2 

-  5.8 

-14.75 

TABLE   1 
PERCENT  OF   ERROR  INyV  AND  )    BY  USING  FIRST  APPROXIMATION 
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<r  : 

=  200 

<r,  ; 

=  500 

M 

f 

453 

2504 

431 

1146 

432 

719 

440 

508 

438 

369 

438 

256 

435 

150 

438 

84 

<r~  ■ 

=  500 

<T7  ■ 

=  200 

-Q-  /A  f  PKILL 

10  453  2504  .419 

20  431  1146  .582 

30  432  719  .652 

40  440  508  .684 

50  438  369  .692 

60  438  256  .699 

70  435  150  .712 

80  438  84  .716 


/*  f  PKILL 


10 

365 

2002 

.227 

20 

315 

769 

.402 

30 

391 

658 

.499 

40 

459 

661 

.569 

50 

472 

538 

.615 

60 

468 

429 

.637 

70 

463 

306 

.642 

80 

471 

155 

.649 

TABLE  2 
PROBABILITY  OF  DESTRUCTION  AS  A  FUNCTION  OF  RE-ENTRY 
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obtained  by  the  use  of  dividers  applied  to  graphs  like  Fig.  5. 
The  results  are  corrected  by  curves  like  Fig.  6. 

As  was  expected  the  method  of  section  III  gives  a  more 
accurate  first  approximation  when  the  ballistic  sigma  is  small. 
Probabilities  of  destruction  can  be  predicted  more  accurately 
for  higher  values  of  6".   In  general,  the  method  of  section  IV 
gives  values  of  AAand    J   that  are  low  and  nearer  to  target  for 
small  -6>"(©-  —  40°)  ;  high  and  further  from  target  for  higher  -Q- . 

This  has  high  intuitive  appeal  since  if  re-entry  is  at  a 
small  angle  we  tend  to  shoot  for  the  point  in  Fig.  5  where  the 
destructive  envelope  intersects  the  ground.   But  if  ballistic 
errors  are  allowed,  we  run  a  high  risk  of  having  the  vehicle 
impact  with  the  earth  before  ever  reaching  target.   Similarly, 
for  steep  angles  of  re-entry  we  tend  to  aim  our  errorless 
vehicle  to  pass  thru  the  highest  point  on  the  envelope  since 
the  vertical  distance  traversed  is  a  maximum.   But  if  we  allow 
errors,  we  run  a  high  risk  of  passing  to  the  left  of  the 
envelope  of  Fig.  5. 

Appendix  C  contains  a  list  of  starting  points  (AA  ,   f.   ) 

f or  ■£*  =5,  6,  7,  ,  90  along  with  associated  destruction 

probabilities.   If  one  chooses  to  use  the  curves  given  in 
Fig.  5,  the  destruction  probabilities  given  may  serve  as  a 
check  on  the  final  probabilities,  the  percent  of  difference 
decreasing  as  (ji    decreases. 
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Outside  of  Appendix  C,  all  probabilities  were  computed 
using  the  polynomial  fit  to  the  curve  since  the  curves  were 
not  believed  to  be  accurate  in  the  first  place  and  we  were 
interested  primarily  in  technique. 

Graphs  of  the  polynomial  used  and  the  actual  curves  are 
found  in  Fig.  7 . 
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FIGURE  7 

PLOT  0?  ACTUAL  DESTRUCT  CURVE  AND  FIT 

(both  scales  in  hundreds  of  feet) 


V 
CONCLUSIONS 

The  optimum  probabilities  of  target  destruction,  and  the 
optimum  values  of  A*\  and  j  appear  as  a  function  of  &  in 
Table  2  for  parameter  (<Tf  <Tt  )   =(200,  500)  and  (500,  200)  ft. 

Notice  that  the  fusing  error  parameter  (P  is  the  more 
important  parameter  in  terms  of  increasing  the  probabilities. 

This  paper  should  be  viewed  as  a  pilot  study.   The  method 
should  be  sharpened  in  the  following  ways: 

1.  The  overpressure  envelope  curves  (e.g.,  Fig  2),  should  be 
established  as  a  function  of  yield  and  mathematical  methods  of 
characterizing  them  (e.g.,  Fig  7),  should  be  improved. 

2.  Sharper  numerical  integration  techniques  should  be  applied. 
The  present  Monte  Carlo  method  has  a  3  standard  deviation 
tolerance  of  +  .015  for  the  probability  of  destruction. 

3.  The  technique  should  be  applied  to  a  large  number  of  yields 
using  realistic  delivery  errors  and  fusing  parameters  (  (Ft  and  (f*  ) 

4.  The  hand  method  for  operational  computations  needs  proper 
development,  particularly  in  the  area  of  sharpening  the 
correction  curves  (e.g.,  Fig  6  and  Table  1). 

For  the  parameter  values  used  in  this  study,  the  probability 
of  destruction  was  not  very  sensitive  to  changes  in  the  values 
of M  and  F  .   This  could  change  with  other  parameters.   For 
example,  a  5  MT  weapon  delivered  from  2000  miles 
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may  well  have  a  delivery  error  of  say  CTj  =  5  miles.   If  50  P.S.I, 
is  required  for  target  destruction,  it  may  be  necessary  to 
determine  AA  and  f~  very  precisely. 

It  was  felt  the  assumption  of  a  straight  line  trajectory 
is  fairly  realistic  and  that  every  moderately  large  deviation 
from  a  straight  line  will  yield  approximately  the  same  results 
since  the  major  dependency  is  simply  in  the  fusing  0"~  which  is 
largely  dependent  on  the  particular  warhead  and  not  on  the 
trajectory.   Since  curved  trajectories  mean  lower  speeds,  the 
fusing  error  ^will  be  smaller.   It  follows  that  the  probability 
of  destruction  will  be  higher.   Thus,  the  use  of  straight  line 
trajectories  will  yield  probabilities  which  may  be  used  as 
lower  bounds. 

The  normalicy  assumption  in  the  fusing  (T  is  felt  to  be 
pretty  good  if  the  warhead  is  fused  by  a  timing  mechanism  or 
a  radar,  but  barometric  devices  are  notoriously  non-normal. 

The  circular  normal  distribution  assumed  for  the  trajectory 
is  probably  not  realistic  but  is  of  little  consequence  since  the 
model  is  broad  enough  to  include  different  cross  range  and  down 
range  sigma's  if  necessary. 

The  definite  range  law,  as  explained  in  the  introduction, 
generally  gives  useful  results  and  we  feel  it  is  quite 
appropriate  for  this  problem.   It  may  further  be  noted  that 
this  problem  would  not  be  solved  were  it  not  for  this 
assumption. 
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APPENDIX  A 
A  brief  description  of  the  integration  technique  employed. 

Again  consider  a  20  P.S.I,  curve  in  3-space.   We  will  first 
enclose  the  surface  in  a  "box",  then  choose,  at  random,  a 
coordinate  inside  this  box.   If  the  coordinate  is  inside  the 
surface  of  revolution  we  record  the  "height"  of  the  function 
at  that  point;  otherwise,  we  record  a  zero.   The  sum  of  these 
observations  divided  by  the  number  of  observations  will  give 
the  "average  height"  of  the  function  over  the  surface.   The 
average  height  multiplied  by  the  volume  of  the  box  will 
approximately  be  equal  to  the  integral  of  the  function  con- 
cerned.  As  the  number  of  trials  gets  large  the  approximation 
will,  of  course,  be  better. 

If  we  consider  each  coordinate  chosen  as  a  Bernoulli 
trial,  the  "worst"  the  variance  could  be  is  (J~  ~    Jl_2"   where 
p   =   Pr   (coordinate  is  under  surface) 
q   =  1-p 
n  =  #  of  trials 
The  standard  deviation  is  (f 
In  our  case  we  chose  n  =  10,000  which  yields  an  upper  bound 
(p  =  q  =  h)    on  tne  deviation  as\/%  :  \    :  io  000  =  -005  or  a 
3  (f  'accuracy1  of  .015 
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Xi  : 

Uniform  (-  c,  c  ) 

X2  : 

Uniform  (-  c,  c  ) 

X3  : 

Uniform  (  0,  b  ) 

We  chose : 
Y\ ,  Y2  ,  Y3=  Uniform  (0,1)  random  variables 

XX  =  2  c(Yi-  0.5) 

X2  =  2  c(Y2-  0.5) 

X3  =   b  Y3 

(  b  =  1000,  c  =  2000  in  our  case) 

If  the  triple  of  random  numbers  (X\ ,  X2,  X3)  drawn  satisfy 

2    2    2 
the  inequality  X1  +  X  <  g  (Xo)the  function  is  evaluated  at 

the  point  (X^ ,  X2 ,  X3) ,   otherwise  a  zero  is  recorded. 

Let 

2    2       2 
h(X±)      =  integrand  evaluated  at  (Xi ,  X2,  X3)  if  X1  +  X2  -   g  (X3) 

.     ^,  otherwise 

we  know  tha  t 

li*      1    r    h  (Xi)  -j>  FfAix)]'  f^W^ 

where  ^X)  =  1  in  our  case.   The  required  integral  is  therefore 

1    * 
approximated  by        2   —    "^* 

be    n    ^_  "h  (Xi) 

USE  OF  PROGRAM 
The  following  procedure  is  recommended  for  finding  M.f and  PKILL: 

1.  Function  X  GOF  must  be  rewritten  to  fit  destruction  curves 
employed. 

2.  Values  of(pand(T)   must  be  changed  in  the  main  program  to 
suit  the  weapon  employed. 
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3.   Since  only  curves  for  a  1KT  weapon  are  generally  available, 
the  problem  must  first  be  solved  for  this  case  using  as  a  first 
1  approximation  the  value  of  (M ,  f  )  found  by  the  method  of 
section  IV   and  the  parameters/^  f   and  (f   must  be  multiplied  by 
W     ((J*  is  not  a  function  of  the  yield  since  it  is  inherent 
in  the  fusing  mechanism) . 

Conditional  probability  of  destruction  given  the  "best" 
trajectory  when 

THETA  =  angle  of  re-entry 

A  =  lower  intercept  of  trajectory  with  destruction  envelope 

ALPHA  =  upper  intercept  of  trajectory  with  destruction  envelope 

Umax  =  vertical  coordinate  of  IBP 

XCOORD  -  down  range  horizontal  coordinate  of  IBP 

J  =  down  range  distance  from  target  the  trajectory 
intercepts  the  ground 
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APPENDIX  B 

PROGRAM  PROB 

TYPE  REAL  MU 

DIMENSION  MU(50),SI(50) 

PRINT  11 
11  FORMAT( 1H1  ) 

PRINT  14 
14    F0RHAT(4Xi 1  HI, 5X.5HMU  ,5X,  5HSI  ,10HEXPOFXl  ,10HEXPOFX3 


1.10HPKILL  ) 

SI  ( 1 )    =    670.    $    MU(1)    =    425.    $1    =    1 

A=.(2.*3.  1415926535)»*1.5 

B=l./A 

152  THETA    *    U0./57. 2957795131 

153  SIGMA    =200.  $    F    =SIGMA*SIGMA       $    G=l./F 
CC=TANF(THETA)  $       C    =    l./CC 

S1GMA1    =    500. 

D  =  SIGMA1»SIGMA1 

E  =  l./O 

DO  777  LL  =  1,70 

N  =  0 

VS1  =  VS2  =  VS3  =  0. 

160  CONTINUE 

Zl     =    RANF(-l)     $    Z2    =    RANF(-2)    $    Z3    =    RANF (-3) 
XI    =    2000.MZ1-.5)    $    X2    =    20  00.*  ( Z2-.5)    $    X3    =    1000. *Z3 
Yl    =    X1*X1    $    Y2    =    X2»X2    $    Y3    =    X3*X3 
1023    W=XG0F(X3) 

ww=w*w 

IF(Y1+Y2.LE.WW)  101,102 
102  V  =  VI  =  V2  =  0.     $  GO  TO  575 
101  CONTINUE 

H  =SIGMA*D 

P=1./H$Q=B*P 

161  RR  =  Xl-SK  I  )+X3*C 
R=RR*RR$S=X2*X2$TT=X3-MU( I)$T=TT»TT 

167  U=EXPF(-.5*UE*(  S  +  R)  )+T«G)  ) 

V  =  Q*U  $  VI  =  Xl«V  $  V2  =  X3*V 
575  CONTINUE 

VS1  =  VS1+V  $  VS2  =  VS2+V1  $  VS3  =  VS3+V2 

N  =  N+l 

IFIN.LT. 10000)  160,500 
500   XN=N 

PKILL  =  40000000C0.*VS1/XN 

El  =  4000000000. »VS2/XN   | 

E3  =  4000000000. *VS3/XN 

PRINT  15, I,MU(I)  ,SI(  I),E1,E3,PKILL 
15  FORMAT ( 1X,I4,F10.5,F10.5,F1D.5,F10.5,F10.5> 

I  =  1  +  1 

MUU)  =  E3/PKILL 

SKI)  =  (El+MUtI  )*PKILL*C)/»KILL 
777  CONTINUE 

END 


FUNCTION  XGOF  (Y) 
1022  Bl=  71.1229667E  CI 
B2=-31.7485590E-C1 
B3=  66.9311997E-C3 
B4=-46.2518018E-C5 
B5=  14.8951096E-C7 
B6=-24.1298172E-10 
B7=  19.0555101E-13 
B8=-58.4969877E-17 

XGOF  =   (((((((  B8*Y)+B7)»Y+B6)»Y+B5)*Y+B4)»Y+B3)*Y+B2)*Y+B1 
END 


APPENDIX  C 


THETA 

I. 
■I 

9 

10 


ii 

12 

13 
14 

15_ 

lo 

17 

18 

19 
20_ 
"21 

22 
23 
24 

25 

2<T 

2? 

28 
29 

3  r 

32 
33 
34 
-35 


N37 
38 
39 
40 
ET 
42 

43 

44 

47 
48 
49 
50 


499.999 
424.907 
189.185 
369.360 

351.793 
329.999 
302T32T 
269.999 
235.552 
201.599 
169.999 
"14*1737^ 

114.999 
89.067 
61.310 

29_l2_9JL 

""•00*t)0u 

-00*000 
-00*000 

-00*000 

QQ.QQO 


V 


ALPHA 

.659.945 
660.099 
653.707 
646.386 

640.607 
636.991 


634.545 
634.304 
633.962 
633.26I 

631.038 
630.126 

629.931 
630.912 

637.916 
644.412 
652.994 

663.553 
TT5T8-59- 
689.581 
704.321 
719.641 
735.102 

764.841 
778.468 
790.967 
80_2.227 
"8127226"' 
821.032 
828.780 
835.665 
841.917 
W777B- 
853.488 
859.259 
865.257 
871.596 
"F7373T8" 
885.401 
892.?48 
900.205 
907.564 


Jmax 

587 
531 

505 

494 

487 

479 

"4o8~ 

454 

438 

420 

401 


1W 
36.9 
356 
346 

235L 

3y- 

333 
334 
336 

■JUL 

3^0 

352 

359 
366 
374 

"3BT 
388 

395 
402 
408 

Tiy 

421 
427 
432 
437 

446 
450 
454 

457 
~%o(T 
462 
463 
465 
466 


J; 

XCOORD 

30.023 
60. 080 
46.210 
(  1.486) 
47.129 
77.037 

"~90."5crar" 

92.987 
91. 381 
91.176 
95.059 

TGZTrzer 

111.540 

117.810 

118.230 
111.320 

9oT4'4U"~" 
84.139 

75M5 

80.023 

102.240 

T3779U-0-* 

166.230 

139.500  - 

566 

_932 

—  >'0  •  0^3 

•  103.780 

.117.540 

.131.330 
;i_44.830_ 

"i~5'/TBY0 
•170.220 
■  182. 160 
■193.820 
•205.430  . 


-6e 


1 i- 


-2177220" 

-229.370 

-241.960 

-254.920 

-263.120 

■=28TT25Cr 

-293.940 

-305.720 

-316.090 

-324.550 


.  h 

■4899 
•3865 
.3178 
.2712 
•  2332 


-1923 
•1738 

-1564 
-1399 

-1096 

-964 
-849 
-753 

-620 

-533 
-562 
-560 
T5ST 
-5?6 
-596 
-620 


ro7o" 

•687 

•704 

.716 

■725 

729" 

730 

•728 

•724 

■719 


-7w 

-709 

•705 

.702 

.700 

~<s?g- 

.700 

.700 

.700 

.700 
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APPENDIX  C  (Con't) 


THETA  PROBABILITY  OF  DESTRUCTION 

SIGMA  «   10     50     100     '  200     500 


1000 


1                  u 

000          .890        .576 

.311 

.127 

.064 

6 

.981        .760 

.443 

.186 

.09^ 

7 

.992        .814 

.492 

.209 

.105 

8 

.99*+       .834 

.511 

.218 

.110 

9 

•■$94       .851 

^530 

.227 

.115 

10 

r^ 

>98 J? 

>99        .9 

75 

XW7    • 

9  AU      .. 

■    \ 1 ?? — 

11 

%s 

04 

.595 

.261 

.132 

.145 

12 

1.000      .932 

.633 

.285 

S 

.95^ 

.681 

.310 

.158 

.969 

.720 

.335 

'    .171 

15 

.Q7Q 

.9«tt 

.357 

— 4j-a3__ 

16 

.986 

.780 

.376 

1,1  ^^^^ 

.194 

17 

.990 

.803 
.824 

.39^ 
.412 

.204 

18 

.993 

'.213 

19 

.996 

.845 

.430 

■  .224 

20 

.QQ7 

.86?    ... 

.452 

.??6 

21 

.999 

.889 

.477 

.250 

.    22 

.999 

.909 

.502 

.265 

u 

1.000 

.923 

.521 

.277 

1.000 

.925 

.524 

.278 

25 

.999 

.90? 

,    .4Q3     . 

.7^0 

26 

.999 

.909 

•     .501 

.265 

27 

.999 

.915 

,510 

.270 

28 

1.000 

.922 

.519 

.275 

29 

.928 

.528 

.281 

30 

/q?4 

.  <^R 

.OR? 

31 

.939 

.5^7 

.292 

32 

.9^- 

.556 

.298 

3 

.948 

.564 

.303 

-  .952 

.571 

.308 

35 

.655 

.  <7fi 

.  119 

36X 

.958 

.583 

.315 

37 

.960 

.588. 

.319 

38 

.962 

.593 

.321 

40 

.963 

.5977 

.324 

.Q65 

^nn 

•  •>oA 

.966 

.603 

.328 

42 

1,   m 

.967 

.607 

.330 

44 
^5 

.968 

.610 

.323 

.969 

.613 

•  335 

.071 

.A17 

•  «a«a-7 

47    . 

.972 

.620 

.339 

.973 
.974 

.624 

.3^2 

48 

.628 

.3^5 

49 
if  A 

■ 

.976 

.632 

.3^7 

50                                                                     y 

f          \ 

'              1 

' 

.977 

.636 

.350 
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APPENDIX  C    (Con't) 

THETA  A                   ALPHA             /Jmax  XCOORD  ^ 

51  00.000  909.663             466  -330  695 

52  911.76?          .  467  -333  689 

53  921.521  •  467  .33^  681 
5^  930.684  467  -331  669 
55  932.232            468  -324  653 

3?  933-000     4&8  -315  633" 

57  933.000  468  -303  611 

58  467  -289  585 

59  467  -273      •  356 

60  467  -256  526 

62  466  -223  466 

63  465  -207  438 

64  465  -194  412 


65  464  -I83  389 

£5  464  -174  371 

67  464  -169      •  356 

68  464  -165  345 

69  i 465  -164  .  338 

70  467  -164  •  333 

71  !  ZJ70  Ti64     '  330" 

72  474  -165  326 

73  ^79  -164  322 

74  484  -160  314 

75  490  -154  301 

75  494-  TT2+4  283 

77  495  -130  259 

78  >.  493  -113  230 

79  484  -937  196 

80  468  -732  160 

Si        !  TO  I53J  iz6 

82  401  -374  986' 

83  350  -269  801 

84  288  -241  756 

85  221  -297  865 

H5  135  1530  rrr 

87  116  -588  142 

88  124  -682  161 

89  221  -538  136 

90  468  -105  165 
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APPENDIX  C  (Con't) 


THETA 


SIGMA 


PROBABILITY  OF  DESTRUCTION 
10    50    100    200    500    1000 


51                                1.000  1.000  1.000 

.978 

.640 

.353 

52 

.979 

.643 

.355 

is 

.979 

.646 

.357 

/ 

.980 

.648 

.358     • 

55 

.649 

.3*59 

56 

.649 

57 

.■ 

- 

58 

* 

V 

59 

60 

61 

1 

62 

63 

l 

64 

* 

65 

00 

' 

67 

- 

\ 

68  ■ 

69 

70 

1 

71 

• 

72 

■»... 

73 

. 

74 

75 

76/.: 

' 

77 

78 

79     • 

80 

a 

1 

82 

* 

; 

83 

84 

85 

/ 

gb 

1 

• 

• 

87      ' 

I 

88 

1 

89 

1 

90              .                           \ 

1      \ 

/-    \ 

' 

V 

p 

\ 

(      ' 

s 

f 

/ 
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